9674 vovN 



NATIONAL ADVISORY COMMITTEE "S 
FOR AERONAUTICS 


TECHNICAL NOTE 3386 


SOME CONSIDERATIONS ON TWO-DIMENSIONAL THIN AIRFOILS 

DEFORMING IN SUPERSONIC FLOW 

By Eugene Migotsky 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 



Washington 
January 1955 


] 



TECH LIBRARY KAFB, NM 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


tech library KAFB, NM 



OOtbMST 


TECHNICAL NOTE 3386 

SOME CONSIDERATIONS ON TWO-DIMENSIONAL THIN AIRFOILS 
DEFORMING IN SUPERSONIC PLOW 
By Eugene Mlgotsky 


SUMMARY 


The aerodynamic characteristics of indicially canibered two- 
dimensional airfoils in supersonic flow are determined theoretically. 

These indicial functions are used to determine the power required to sus- 
tain a general time-varying chordwise deformation. The harmonically 
oscillating parabolic mode is considered in detail and stability bounda- 
ries are presented for this case. Also, the thickness distribution of a 
beam having a parabolic fundamental bending mode, in vacuo, is determined. 


INTRODUCTION 


Until recently, little consideration had been given to the aero- 
elastic phenomena associated with chordwise deformations of wings. Thus, 
while spanwiee elastic modes were considered in most analyses of aero- 
elastic phenomena, it was generally assumed that the displacement of a 
spanwlse station was restricted to translations and rotations. In recent 
years more attention has been given to the problems associated with chord- 
wise deformations. One of these problems is the chordwise bending of a 
two-dimensional airfoil in a supersonic stream. 

Blot (ref. 1) appears to have been the first to consider the static 
Instability of a two-dimensional airfoil in a supersonic stream. Miles, 
in reference 2, investigated the static and dynamic stability of a can- 
tilevered (chordwise) airfoil section in supersonic flow by including 
low-frequency aei*odynamlc forces. In the present paper we consider the 
dynamic behavior of a two-dimensional airfoil in a supersonic stream, 
without the restriction of low frequency in the determination of the air 
forces. 

The aerodynamic characteristics of airfoils which are indicially 
cambered in a relatively arbitrary manner are deteimined, by linear 
theory, from known, or easily obtained, characteristics of an indicially 
plunged airfoil.^ These indicial characteristics are then xised to obt ain 

^In work performed independently and concurrently, Lomax, Fuller, 
and Sluder (ref. 3) have obtained similar results for the rectangular 
wing, but by a different method. 
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the rate of energy input (povrer) required to sustain a general time- 
varying chordvise deformation. 

The particular case of the harmonically oscillating parabolic arc 
vas chosefn for detailed consideration. In addition, the thickness distri 
bution of a beam having a quadratic normal bending mode, in vacuo, was 
determined. 


SYMBOLS 


ao 

Ao 

Ah, Ah 

c 

Co 

Cm 


half-amplitude of oscillation of leading edge 
section mean-power coefficient (See eq. ( 2 l)*) 
functions defined by equation (15) 
chord of wing 

speed of spund in free stream 

section lift 


section lift coefficient. 


q.oC 


section pitching -moment coefficient, ^ 

measured about leading edge 


X section moment 





Cma 


cmq 


^,n 

E 

El 


indicial section lift coefficient due to angle-of-attack 
change only 

indicial section lift coefficient due to pitching on a 
wing rotating about its leading edge 

indicial section pitching-moment coefficient due to angle- 
of-attack change only, measured about the leading edge 
and considered positive when the trailing edge is forced 
down 

indicial section pitching-moment coefficient due to pitch- 
ing on a wing rotating about its leading edge, measured 
about the leading edge and considered positive when the 
trailing edge is forced down 

generalized indicial force coefficient (See eq. ( 6 ).) 

Young* s modulus 

E 

1 - 
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F,G 

G-BQ^Gs3^,Ggg 

h 


Ik) 

I 

Jn 

k 

M 

P 

Ap 

P 


Pn 


Pli^Pla 


q. 


<lo 

t 

ta 

'*'0 


Sctcwartz functloa, + ifoj 



functions defined ty equation (18) 
functions defined "by eq,uation (a8) 

functions defined by equation (A 9 ) 


function giving timewise dependence of vertical position 
of wing 

half-amplitude of oscillation of midchord point measured 
relative to leading edge of wing (sketch (c)) 

moment of inertia of cross-sectional area of beam with 
respect to neutral axis 

Bessel function of the first kind of order n 


reduced frequency, 


(jJC 

2Vo 


flight Mach number 


static pressure 


discontinuity in pressure across z = 0 plane 

loading coefficient, ^ 

qo 


loading coefficient correspondlng^to a normal velocity 
distribution equal to - 


loading coefficients corresponding to unit steps in h(t) 
and h(t) , respectively 

0c 

Vo 

1 2 

dynamic pressure, -g PqVq 
time 

time to reach steady state 

Cpt 

c 



k 
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c^t 

* 

T 

maximum thlcloaess of airfoil section 

A 

Vo 

flight velocity 


Wu 

pertiirbation velocity component in z direction on upper 
surface of z = 0 plane 

— 

w 

energy required to sustain the motion 


x,z 

Cartesian coordinates 


a 

angle of attack, radians 


r 

Gamma function 

— 

? 

function giving chordwise dependence of vertical position 
of wing section 


n 

thickness of beam 


e 

angle of pitch, radians 

A 


mass per mit length 


V 

Poisson's ratio 

V 

1 

distance from leading edge of wing section 


Po 

density of mdistxirbed air 


Pvr * 

density of wing material 


T 

variable of integration 


(P 

Vot 

chord lengths traveled, 


tPa 

chord lengths traveled to reach steady state 


$ 

perturbation velocity potential 


(0 

frequency of oscillation, radians/sec 


2o 

2kM® 

reduced frequency parameter, 


(*),(') 

differentiation with respect to t and respectively 
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AERODYNAMIC CHARACTERISTICS OF INDICIALLY CAMBERED AIRFOILS 


The probloa considered In this section is the determination of the 
Indicial functions (pressure, lift, and moment coefficients) for a thin, 
two-dimensional airfoil that starts from rest and suddenly moves with a 
constant, supersonic, forward velocity and simultaneously attains a chord 
wise variation of normal velocity. The basic, linearized, partial differ 
ential equation for this problem, in terms of the perturbation potential 
$, may be written 


®xx ®zz ^t^ti 


( 1 ) 


when the fluid at infinity is at rest 
with respect to the xz coordinate 
system. The boundary condition to be 
satisfied is that the normal velocity 
Wu be a specified function of dis- 
tance from the leading edge of the 
airfoil (x + Mtx) in the region traced 
out by the airfoil in the xzti space 
(sketch (a)). The loading coefficient 
P, in terms of the perturbation poten- 
tial, is then given by the relation 


JlJl 

VoM ^ti 


( 2 ) 


Loading Coefficient 



The solution to this boundary- 
value problem has been obtained for 
the plunging airfoil, that is, 

Wu = constant (see ref. 4). At super- 
sonic speeds, this solution applies as 
well to the case shown in sketch (b) 
where Wu = 0 for 0 < x < - Mti + 
Thus, if Po(x,ti) is the loeidlng coef- 
ficient for Wu = 1 in the shaded 
region of sketch (a), and if P(x,|,ti) 
is the loading coefficient for w^ as 
in sketch (b) , then 


Sketch (a) 
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=• Po(x - |,ti) for - Mti + 1 < X < c 
■ = 0 for 0 < X - Mti + I 


Siaoe the theory is linear, the loading coefficient for an airfoil with 
arbitrary Wu_(|) can be obtained by superposition of elementary solutions 
of the form given in equation ( 3 )* Note that in this superposition it 
is necessary that w be a function of | only. “Hie superposition gives 

px+Mti 

P(x,ti) = Po(x,ti) + J w»(|) Po(x - |,ti) d| ( 1 +a) 

which, upon introducing q = x - may be written in an alternate form 
of the Duhamel integral 


P(x,ti) = Po(x,ti) w-uL M+ + r w'(x - cf) Po(d^ti) dff (Ub) 

-'-Mti 


Generalized Indicial Pprce Coefficients 


In characterizing the aerodynamic properties of an airfoil that is 
only pitching and plunging, it is generally sufficient to give only the 
usual lift and moment derivatives (i.e., chlq, etc.). For an air- 

foil undergoing more complex chordwise deformations, it is formd that 
more information about the distribution of loading on the airfoil is 
needed. This additional information is conveniently represented in the 
form of nondimensional generalized force coefficients which are defined 
as follows: Let the normal velocity of the airfoil be given in the form 

n n . 

and let Pn be the loading coefficient corresponding to a normal veloc- 
ity distribution equal to (. ? an = 1, Sm = 0 

for m n. We define the generalized force coefficient as 


Cm,n 



dx 


( 6 ) 
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T 




If 


For the special cases of m and n equal to 0 and 1 ^ equation (6) reduces 
to 


■' 0,0 


-0,1 


= -c. 


^CL 


= -c 




Ci^o - 

= Cniq 


(7) 


Thus, the coefficient Cm^n ^7 considered to he a generalization of 
the usual lift and moment derivatives . 

Inasmuch as the loading corresponding to an arbitrary normal-velocity 
distribution can be expressed as a function of the uniform-downvash load- 
ing, it follovs that the force coefficients for arbitrary n can be 
obtained from those corresponding to n = 0 * To this end "vre substitute 
equation ( 4 ) into equation (6) to obtain, for n > 0 (so that 


_ Vo + Mti' 


Vo 

cm,n = -5- / V* 

° J _Mti V 


tTax r £( 

/ ^ -Mti 


(x - Po(cr,ti) da 

on 


Interchanging the order of inte^ation, ve have 
p -Mt i+c p -Mt i+c 


p — Mt i+c P -Mt i+c 

/ Po(cr,ti) dcr / (x + (x + dx 

'•'a 


Denoting the inner integral by l(ff) and expanding (x - cr) by the 
binomial expansion ve obtain 


1(b) = c 


n+m 


/n - IN (-1)^ /g + MtA ^ 
lNr^m + n- r\ c J 

r=o 


r=o 
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The second integration then gives the desired recxirsion relation 
for n > 0 





(-1)“ Cm+n,o ^ ^ 


- 1 )^ 
r + 1 


vhich relates tlie. coefficients for arbitrary n to those for n = 0*^ 
The coefficients Cm o indicial plunging are given in the Appendix 

for m between 0 ana and^ in integral form^ for arbitrary m. 


POWER REQUIRED TO SUSTAIR A CHOEDWISE DEFORMATION 
General Considerations 


The dynamic stability of a linear oscillating system is determined 
by the rate of energy flux into the system# For an airfoil deforming in » 

an air stream^ this energy arises from the aerodynamic forces acting upon 

the airfoil; if power is req,ulred to sustain the motion against the lift 

forces, the motion is stable, conversely, if power must be extracted from 
the system in order to sustain the motion, the motion is unstable. 

Let us assume that the position of the mean line z, at any instant — 
of time is given by (in an axis system attached to the airfoil) 


z = h(t) Hi) (9) 


Then, the power required to sustain this motion against the lift forces 

W = - J Ap(|,t) || d| 


®An equivalent recurrence relation was obtained In reference 3 foi* 
the more difficult case of the rectansrular wing. 
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33 iay be ■written 


W = - h(t) r Ap(|,t) ^(1) d| 


( 10 ) 


The normal velocity distribution across the airfoil, which deter- 
mines the loading on the airfoil, is given by 


8z . 8z 


^u - 


( 11 ) 


or 

Wu = Vo^’(l) b(t) + ^(1) h(t) (12) 


This downwash distribution may be considered as the sum of two contribu- 
tions j the first arising from the instantaneous position of the airfoil 
and the latter resulting from the instantaneous normal velocity of the 
airfoil relative to an observer moving with the airfoil. Now the load- 
ings resulting from each of the terms in equation (l2) may be obtained, 
by the Duhamel integral, from the loadings corresponding to unit steps 
in h(t) and £(t) . Thus, the loading coefficient for a motion as des- 
cribed by equation (9) ts 

P(|,t) = Ph(^>'^) ~ ^ 


where Ph is the loading coefficient corresponding to a unit step in 
h(t) , that is, Wu = Vq^'U), and P£ is the loading corresponding to a 
unit step in ]a(t) , that is, W'u^ = ^(^). 

Substituting equation (13) into equation (lO), and interchanging the 
order of Integration, we obtain the following expression for the power 
required to sustain the motion against the lift forces 


W = - 



t) Ah(T) dT + h(o) Ah(t) + 


c 




t) A]^/t) 


dr + :^ h(o) 



(IJ^) 
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where 


Ah(t) 

Ah(t) 



Ph(i>t) ^(1) 




= ^ r°Ph(i,t) ^(1) 


(15) 


The indiclal loading coefficients Pj^ and P£ can he obtained from 
equation (if-) and the problem of obtaining the power required to sustain 
a given motion is reduced to a series of quadratures. It shoiild be noted 
that if the function ^(|) is a polynomial, the functions A^ and A£ 
can then be obtained in terms of the generalized force coefficients Cm n* 
This will be Illustrated later where is a quadratic function. 


Application to Sinusoidal Oscillations 


We assume now that h(t) is given by 


h(t) = ho sin ut 


(16) 


Substituting into equation (lif-) gives 


W 


qpOJ^ho^ 

2 


a 


All '^.T + 


IOC 

Vo 


I 


' 2c 

Afi Bin wt dr + ^ Afi cos 0 )t + 
Vo 


a 


it p t \ 

sin 2cot ( / Aji sin cot dr - j A;^ cos cot dr^ + 


cos 2cot 


(X' 


(OC 


Ah cos COT dr + ^ / Ah sin cot dT 


Ah sin COT dT^ 


( 17 ) 


Introducing the functions F(t) and G(t) such that for 0 t ^ t: 


Ah(T) = Ah(t) + F(t) 
Ah(T) = Ah(t) + G(t) 


( 18 ) 
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and substituting into equation (17) gives, after integrating and simpli- 
fying. 


= 


qobo^^ 

2 



COS COT dr 



sin COT dT + 


cos 


sin Scot 


Scot F(t) cos cot dT + ^ 

t) sin OJT dr - cos 0)T (19) 



Since the functions F(t) and G(t) vanish for t = ta, vhere ta is the 
time to reach steady state, the integrals in equation (19) become con- 
stants. Thus, the power required to sustain the oscillations, after the 
transients die out, reduces to the form 


W = 


qoho^co^ 

2 


(a + A cos 2 ut + B sin 2 cat) 


( 20 ) 


Upon substitution of the nondlmensional frequency 
lengths traveled cp = equation (SO) becomes 


coc 

2 Vo 


and the chord 


^ ^ 2ho^k^Voqo 


(Ao + Ao cos 4 kcp + Bq sin it-lop) 


(21) 


vhere 


P 'Pa p Pa 'I 

Ao ~ J ^(p) cos 2kcp dip - A^(pg^) - U(p) sin 2kcp dq) 


Bo ■ 


Aht-pa) . 

2 k 


o 
9 a 


pcpa , ^ P 9 a 

/ F(cp) sin 2 kcp dcp + k / G(cp) cos 2 kcp dtp 
*^o '-'o 


> ( 22 ) 


and 


F(cp) = Ah(cp) - Ajj(cpa) 
G(cp) = Ah(cp) - A£(cpa) 


Thus W is a function which oscillates about a mean value with twice the 
frequency of the oscillation of the motion. For considerations of dynamic 
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stability, the constant Aq is the significant parameter. A positive 
value of Ao indicates that^ on the average , there must be a steady influx 
of energy into the system in order to sustain the motion, that is, the 
system is stable, conversely, a negative value of A© means instability. 


QUADRATIC CAMBER DEFORMATIONS 


We shall consider, in this section, the application of the preced- 
ing analysis to the case of an airfoil which is assumed to oscillate in 
such a manner that, at any Instant of time, ‘its mean line is a simple 
quadratic function. A question which naturally arises is: Does this 

assumed quadratic function correspond to a nozmial vibratory mode of a 
beam whose thickness distribution is a reasonable approximation to those 
used in airfoils? In order to help answer this question, we shall first 
consider the problem of determining the thickness distribution of a beam 
which has a prescribed vibratory mode shape;® in particular, we shall pre- 
scribe the assumed quadratic function. 


Detezmlnation of Beam Having a Quadratic Normal Bending Mode 


The differential’ equation for the transverse vibrations of a simple 
beam of nonuniform thickness is given by 








(23) 


The boimdary conditions for a free-free beam are given by the vanishing 
of the shear and bending moment at the ends of the beam 

(eiI Ip) = ^ (eiI IJi) =0 at X = o and X = c (2it-) 

Specifying the mode shape to be 


= h^ 


sin cot 




c| 



(25) 


g 

Mathematically, this problem is the inverse of the usual vibrating- 
beam problem in which the beam thickness distribution is specified, 
together with certain boundary conditions, and the normal modes are to be 
determined. Here, in the inverse problem, a mode shape is specified and 
the thickness distribution is to be obtained. 
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and substituting into equation (23) we obtain 


8 eA 

d| 


2 = 


^ (cl - 1^) - ^ 




I 


(26) 


Since, for a solid beam of rectangular cross section with thickness i] 
and width b 


I = and [j. = pwbTi 

equation (26) may be written in the form 

\<i|/ d?^ \ 2 Ei y \ c c^ ho/ 


Similarly, the boundary conditions ( 24 ) reduce to 


Tl(o) = q(c) = 0 


(27) 


( 28 ) 


Assuming a solution of the form 


q = bo + bxl + bsl^ (29) 

substituting into equations (27) and (28)^ and solving for the coeffi- 
cients bo^ bi^ and bg, we find that 

ap _ 4 

bp 5 
bp = 0 

bi = c(x)*yM./5Ei 

ba = - 


or that 


Tj = (J yp^/5Ei (c| - 1^) 


(30) 
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Upon siibstitutlon of the maximum thickness T, the thickness distribution 
may also be written as 


T 1 = -^ (C| - 1 ^) (31) 

From equations (30) and ( 31 )^ it follows that the frequency of the assumed 
mode shape (eq. (25)) must be 

“ " c (c) '/ 5 Ei/pw (32) 

It can be shown further, by applying the Stodola iteration method (see 
"ref. 5, p. 313) that the normal bending mode given by equation (25) 

With 7-^ = ■=■ is the fundamental mode. Thus, we have shown that a solid 

1^0 5 " 

airfoil with a parabolic thickness distribution has, for its fundamental 
chordwise bending mode, one of the assumed quadratic functions. 


Average Power Required to Sustain the 
Quadratic Camber Oscillation 


For 
arc (see 


the particular case of the harmonically oscillating parabolic 
sketch (c)) we have 


h(t) = h^ sin o)t 


^( 1 ) 





(33) 


The indicial loading coefficient 
corresponding to a downwash 
variation 

= U'(i) “ ''o [- (B) ■"©]' 

may be written as 

+ (3W 



Sketch (c) 
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Simileirly, the indlcial loading coefficient corresponding to a down- 

wash distribution w^ = ^(|) = -hU^/c^) + k{^/c) - (ao/ho), becomes 


Ph = -^2 + S (35) 


Substituting equations (33) -(35) iiito equation (l5)j comparing the terms 
obtained with equation ( 6 )^ and using recurrence . relation ( 8 ), we obtain 

3’ *^0,0 ^1^0 " '3 *^3,01) ^ ^ (*^0,0 “ ^^1,0) (36) 

and 

■^h ~ Cx^o - ^ ^ ( - 2cx,o + 2 c2^o) + “^0,0 

(37) 

Substituting the values of Cm q from the Appendix, we obtain Aj^ and 
Ab, as functions of Mach number M and nondimensional time tQj these 
functions are presented in equations (A7) and (a 8 ) of the Appendix, and 
are plotted in figure 1 against cp, the chord lengths traveled, for sev- 
eral Mach ntmbers and two values of ao/ho. 

The mean power coefficient Aq Is obtained by substituting the 
expressions for A-^ and A^ from the Appendix into equations (22) and 
( 23 ) and integrating. The method of integration is very briefly outlined 
in the Appendix and the final results for the Integrals appearing in 
equation (22) are given in equations (A9) and (AIO) of the Appendix. 

The functions Fco^^Cx^^so^tJsx^*^S 2 which occur in equations (A9) and 

(AIO) are plotted in figure 2 against w for several Mach numbers* The 
mean pover coefficient Aqj as a function of ai^ is presented in figure 3 
for several Mach numbers and tvro values of ao/ho- It will be noted frcmi 
figure 3 that regions of instability (Ao < O) exist iat the lower super- 
sonic Mach numbers. 

Whenever regions of both stability and instability exists it is 
usually convenient to delineate these regions by stability boundaries 
which separate the two regions. To this end we note that Aq uiay be 
written in the foim 
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Ao = , + Fco(M^C) + Gso(M, 5 o) 1 + 

|_15 n/m^ - 1 J 


(i 


ao > 

ho. 


16 



3 n/m^ - 1 


+ Fc2^(M,w) + Gs^(M,cS) 


isJW 


== + Gs2(M^“) 


(38) 


that is , a quadratic function of the parameter ao/ho which determines 
the location of the fixed points of the oscillation. Thus, for given M 
and S, at most two roots exist for Aq = 0. For the limiting case of low 
frequency, a simple expression may he obtained for the values of ao/ho 
corresponding to the stability boundaryj this solution is 





+ 12M® - 




(39) 


The stability boiindaries for the harmonically oscillating parabolic arc 
are presented in . figure with ao/ho plotted against £0 for three Mach 

numbers at which the boundaries exist. These boundaries are presented in 
a different form, obtained by cross-plotting the curves of figure 4, in 
figure 5 with ao/ho plotted against Mach number for several values of 
the reduced frequency k. 

It will be noted that the motion is always stable for Mach numbers 
greater.- than I.65 regardless of the location of the nodal points and of 
the frequency of oscillation. In addition, the stability boundaries also 
disappear for all values of Mach number and ao/hq when k is greater 
than 0.65. These limits on Mach number and frequency, when applied to a 
solid steel airfoil with quadratic thickness distribution and with the 
assmption that the mode shape is not affected by the air forces, show 
that the motion can be unstable only if the thickness ratio is less than 
1.5 percent. It should be noted that this limit on thickness ratio is 
valid only to cases where, a one-mode analysis applies. It has been found 
(ref. 2) that a two-mode approximation yields another stability boundary 
which is associated. with the coalescing oflthe frequencies of the two 
modes when the ratio of aerodynamic to elastic forces becomes large. In 
such cases, a more complete analysis would be required. 
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IT 


CONCLUSIONS 


A recTirrerLce relation for the generalized force coefficients for 
Indicially cambered two-dimensional airfoils in a supersonic flow has 
been determined theoretically from the response of an Indicially plunged 
airfoil. The rate of energy input required to sustain a general time- 
dependent chordwise deformation has been obtained in terms of these Indi- 
clal functions. The particular case of a parabolic chordwise bending 
mode was considered in detail, and stability boundaries were obtained 
which show that this mode is mstable only for Mach numbers less than 
1.65 and reduced frequencies less than 0.65* In addition, it has been 
shown that a beam with parabolic thickness distribution has a quadratic 
fundamental chordwise bending mode in vacuo. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 11, 195^ 
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APPENDIX 

GENERALIZED INDICIA! FORCE COEFFICIENTS 


The generalized Indicial force coefficients for a plunging two- 
dimensional airfoil traveling at supersonic speeds may he obtained by 
Integrating eq,uation (6) for u = 0. The results of the Integration are 
given in the following expressions. 


0 ^ to 


< 1 
“ M + 1 


u 


(Ala) 


Ci,o = 



(A2a) 


c 


2,0 


4 , 4to^ 
3M 3 


(A3a) 



(A4a) 


c 


4,0 - 



(A5a) 




4 

(m + 1)M 


r. ni+i pTT jjj 

— — 2 — _ / cos u(M - cos u) du 

ir(m + 1)M Jo 


1 

M + 1 


^ tf 


< 1 
“ M - 1 


(A6a) 


4 

7T 


— arccos 
M 


Mtr 


to 


- 1 


arccos(tQ + M - t^M^) + 


^ sfto® - (1 - Mto)® 


(Alb) 


Cl,0 = - 


V 


^ 1 




arccos 


Mto - 1 


n/P" 


arccos(to + M - toM^) + 



7to^ - (1 - Mto) ^ 


(A2b) 
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^2^0 = 


=3, 


Cm,o 



A (A. 

37T |\M 


Mto - 1 1 

n/m2 - 1 


(|j + f)v] ,/t„= - (1 - 


i + Si2 + (" J. + «'l t 2 

3 M 3 M ^ V 3 M ^ 3 y 


arccos (to + M - toM^) + 

(A3t) 


TT \ 2 y ° 


Mto - 1 

axccos — 4 : + 

■Co 


\/M^ - 1 


arccos (to + M - toM^) + 


■ytc" - (1 - «to)=} 


(A^b) 




arccos 


n/m^ - 1 


arccos (to + M - tc^ 4 ^) + 


Mto - 1 , 
to 


5M 5 


Aj- ^^+3 /^ m£N ^ 3 f _8_ 0^ 

Vf ^ V30 5/ VI5M 30 

to^ Vto^ - (1 - Mto)^| 


5 


Lm 


_4 fl 

+ l)7T 


arccos ~ - + ^ arccos (to + M 


n/M^ - L 


/Mto - A 
, m+x p arccos — =‘r ) 


m 


(M - cos u) du 


(A5b) 

toM^) + 

(A6b) 


- 1 


v/M^ - 1 


(Ale) 


Cl,0 - 


n/M2 - 1 


(A2c) 
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^2,0 “ 


(A3c) 

3 n/M^ - 1 

cs,o = 

1 

(a4c) 

v/m^ - 1 

C4,0 = 

4 

5 n/M^ - 1 

(a5c) 

0 

4 

(a6c) 


(m + 1) - 1 


Fijnctlons Associated With Qiiadratic Camber Deformations 


The indicial coefficients A^. and A£,, for the quadratic camber 
deformation, may be obtained by substituting equations (A 1 ) through (A 5 ) 
into equations (36) and (37) 'to obtain: 


0 <t. < 


o - M + 1 


Av = 16 


¥ - (® - “) - II (- ^ 


A^ = 16 


15M 3M 




(A 8 a) 




M + 1 


M - 1 


An = f (® + + H (- %^)] 


Mto - 1 . 

arccos — ^ + 

to 


~ + 3 )to^ + (13M + ai®)to® + 


(A 7 b) 


(S) - i } 7to^ - (1 - Mto)^ 
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Ah 


fl 1)M2\ 5 

ao^ 

(,5 15 j ° _ 

w ^ 


' TT ■ 15M 3M 3 V5 _ 

^2tof] . /a^Y itM 

3 J vW V M/J 


8to- 

M 


- Mto - 1 . 

arccos — + 

to 


1 r - ^ ~ h ^ 

« 3W WJ 


arccos(to + M - toM^) + 


3“o (f * S *0 ■ - ("5 - 

(^ + 83M + a®) to‘] + (^) [i + ^ - (p + + 

Is) (■ I)} '/to“ - (1 - Mto)“ 


(A8b) 


to ^ 


M 


Ah = 0 


Afi = 



r. ^ H 

H — r~ 

- k (- 

LoVI 

1 

L 15 ' 

3 Vii 

oj Vi 

lo/ J 


(ATc) 

(a8c) 


The integrals appearing in equation (22) for the mean power coeffi- 
cient may te evaluated by integrating by parts, making the transformation 
M 

CP = s u + — = , and reducing the many integrals to the Poisson 

^ - 1 m2 - 1" ® _ 

integral representation of the Bessel functions (ref. 6) 


Jn(z) = 


2 'I 


n/F r 


(“^1) 




cos zu)(l 


u2)^ 2 


By means of successive application of the Bessel function recurrence 
relation 


Jn_j^(z) + Jn+i(z) = — Jn(z) 
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the sine and cosine integrals in equation ( 22 ) can be reduced, after 
considerable manipulation, to the forms 


J4 


pM-i r / \ 

/ F(cp)cos 2kcp dtp = - Fco(M,C) + \J^J Fci(M,w) 

Jq l \ / j 


(A9) 


where 


+ 1 ) 

+ TTT I + 


,r . 1 ^ , /|Y,^ ,r 

L fM2 J VV L 




JqI )cos u| 


- ICM^ 


“ 

+ Jol^lsin Gi 

^\My I /■«2 ^^7/2/'^■.^3 


° [ ( M ^ - ' 1)V^ T 


(m 2 - l)7/2(n)3_ 

-ai®(ai 2 + 13 ) 1 \ 

^^(iFJj 


Fci(M,G) = 16 |fop(M,G)J^- 

J 3 . (|) sin G [ 


- M*= 


(m 2 . 1 )^/^(G)^ 
M® 


(m2 - 1)^/2(£3)2 


:] + Jo (m) <=os 0 _ j^J5/2(0)=J 

] 


and 


I 


M 

M-i 


G(cp)sin 2kcp dtp = 


I j^Gsp(M,G) + (g) Gsi(M.a) + (^g^\g(M,G) 

(AlO) 



I 


vliere 






on / 1 „ 2M^ 12M^ + 16M*^ 1 ^ ^ TaF ^ 

32 -| - ^_ + fOi^C ^^)[ " ,\i=;y!?/-\a “ /,,a «^a/a/•^^ 4 .l ^ ^ ^ 'I " n«#a 


i i;> ^/M^ - 1 


L ^W“ - 1 ;*^' ^jxi“ - 


xpri- 


^ 1 + Jo (g') COB C r j. y — 

3(M= - l) 3 / 2 (w)J W L15 ^/M* - 1 


m 


a 


ISM"* + I6M® 


(h) “ [■ 


15(M® - l)3/'®(0) 15 (M' 


15 (m 2 - 1) 3/2(0) a (M2 - l)®/2(0) 

fflrSv-<%-] • ■■ (j) - ■ [■ ” 


S)'*_ 


15 (m 2 - l)2/2(0) 


1 _ _Q / -<a \ 


15(M2 - 1) 7/2(0) 3, 


J 


y \ 

+ Ji fr 1 sin 


5[- 


^ y _ ^r» _ \ 


k--Ck/A-^ j-v O 


1 - aj ^ gun MS* + baft"- + xo; 

15M*^M2 - 1 15(M2 - l)5/2(o)a 15 (m 2 - l)®/2(0)^ 


■Ji 


. /l^p -1 *1 


^ is'/lFTT ’*' - l)5/2(0)2j S^“'“^1(M2 - 1)3/2(CJ) ^ 3M2 J ■*■ 

w [- ^ /:.F— T “ TTS ivs/a/r. Tal + (f) “ [" o/vre :^n\a/^l + 


Jo IS) cos 
vy 


L J V - X (M" - X/ ' \w; J 


o^urC 

j- - j_/ 


Ji (m) “= ° [3(„2 ”i)=/J ■" (C^ ° [■ sa-j^-rj ■ 3(11=*- 1)=*(0) J } 


1. r 1 „ /„ -\r . jm 2 - i((7iV 


- M ./■ 1 N , / 3 \ . 1 \l 

“Vjp=^ -"W""* 


ro 

CjJ 
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In the limit as frequency goes to zero these integrals become 


lim / 
k— o '^o 


pM-i 
/ F(cp) 


COB 2kcp d£p 


r 

hoj L" 


15(M2 - 1)3/2 \hoy L3(m2 - 1)3/2 


lim k 

k — >0 


pM-l 
/ G(cp) 


sin 2kcp dcp 
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(c) Gs^ and Gg^ 
Figure 2 .- Continued. 
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Figure 2.- Concluded. 
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